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Abstract 



We study the pair production of scalar top quarks ti in e'^e collisions with the 
subsequent decay into 5-quarks and charginos, ii — )■ bxi ■ We simulate this process 
using PYTHIA6.4 for beam energies 2Ebeam = Vs = 350, 400, 500, 800, 1000 GeV. 
/\ • Proposing a set of criteria we obtain a good separation of the signal stop events 

^ . from top quark pair production which is the main background. The number of stop 

production events obtained with the proposed cuts for different energies is calculated 
for an integrated luminosity of 1000 fb~^. We propose a method to reconstruct the 
mass of the top squark, provided the mass of the lightest neutralino is known, and 
estimate the error of the mass determination for the case ^/s = 500 GeV. 

1 Introduction. 

The scalar top quark, the bosonic partner of the top quark, has attracted much 
attention as it is expected to be the lightest colored supersymmetric (SUSY) p] particle. 
tL and tji, the supersymmetric partners of the left-handed and right-handed top quark, 
mix and the resulting two mass eigenstates ti and t2, can have a large mass splitting. It 
is even possible that the lighter eigenstate ti could be lighter than the top quark itself [2], 

Searches for top squarks were performed at LEP and Tevatron and will continue at 
LHC andlLC |1], [5]. 



This Note is the continuation of our previous Note where we have considered stop pair 
production in photon-photon colhsions [6]. In the following we study the reaction 



+ e -> ti + ti 



(1) 



Among the possible ti-decay channels within the MSSM (see [7] for details), we focus 
on the decay ti — )■ bxf followed by the two-body chargino decay xf — ?■ Xi^^^? where one 
of the W's decays hadronically, W — )■ qiQj, and the other one leptonically, W — > yuz/^ \j 
[8]. The final state of this signal process, shown in the left diagram of Fig.l, contains two 
b- quarks, two quarks (originating from the decay of one W boson), a hard muon plus a 
neutrino (from the decay of the other W) and two neutralinos: 



e e~ 



^ iih -^ hhxtXi -^ hhW+W-xlxl ^ bbMj^iP^x'iXi- 



(2) 



The main background process is top quark pair production with the subsequent decay 
t — i- bW^ (for W's we use the same decay channels as in the stop case): 



e~^e -^ tt ^ bbW^W -» bbqiqjiiv^. 



(3) 



The only difference between the final states of stop and top production (shown in the 
right diagram of Fig.l) is that in stop pair production there are two neutralinos which are 
undetectable. Thus, both processes have the same signature: two 6-jets, two jets from W 
decay and a muon. In the following we show that the physical variables constructed of 
jets combinations may allow to reconstruct the scalar top quark mass. 

In the present paper we consider only top pair production as background. We analyse 






a) r M,^ ^^ 

Figure 1: Left is the stop signal event diagram, Right is the top background diagram. 

the processes (2) and (3) with the help of Monte Carlo samples of the corresponding 
events generated by the two programs PYTHIA6.4 [12J and CIRCEl [13]. The program 
CIRCEl is used for the parametrization of the beam spectra involved in the processes (2) 
and (3) to account for the effects of beamstrahlung. The energy of the beams is chosen as 
'2Ebeam = ^ = 350, 400, 500, 800, 1000 GeV. Most plots are shown for ^ = 500 GeV 
and a stop mass Mj = 167.9 GeV. 



The process e+e — > f i + i i with the subsequent decay ti — > cxI was considered in [S] -|11]. 



In Section 2 we give the set of MSSM parameters used in our study. 

In Section 3 we discuss some general characteristics of the signal process e~^e~ — )■ titi 
and the main background e~^e~ — t- ti. First we show how the beam energies are affected by 
beamstralung and other beam interaction effects which are simulated with CIRCEl. The 
subsections include kinematical distributions (obtained without imposing of any cuts) for 
the produced stop quarks, for the jets originating from W boson decay and for 6-jets. We 
compare them in detail with those of top pair production. 

In Section 4 we demonstrate how to discriminate between the signal muons produced 
in W boson decays and those stemming from hadron decays in the same events. 

In Section 5 we show the distributions of the global variables missing energy, total 
visible energy, the scalar sum of the transverse momenta of all visible particles in the 
event and the invariant masses of the final-state hadronic jets plus the signal muon. We 
show that they are good tools for separating the signal from the top background. 

In Section 6 we introduce further two global variables, the invariant mass of the final- 
state hadronic jets and the missing mass and demonstrate that they are very useful for the 
separation of signal stop events from the background top events. We propose three cuts 
which provide a good signal-to-background ratio (S/B). The impact of the proposed cuts 
on the values of the cross sections of stop and top pair productions are shown together 
with the values of S/B ratios for different values of ^/s. 

Section 7 is devoted to the mass reconstruction of the scalar top quark based on the 
distribution of the invariant mass of one 6-jet and the other two non — 6-jets (from W 
decay), provided that the neutralino mass is known. 

In Section 8 we show the distributions of the invariant variables described in Section 
7 for a stop mass Mj^ = 200 GeV. 

Section 9 contains some conclusions. 

2 MSSM parameters and cross section. 

The scalar top quark system is described by the mass matrix (in the ii — iji basis) 

i, M 

f M? M? ^ 



Mf Mf 

tRL tRR 



(4) 



with 



^L = ^1 + (2 " -sin'ew)cos2(5Ml + M^, (5) 

M^ = Ml + -sin^Qwcos2l3Ml + Mf, (6) 

RR 2 

^L = (^D* = ^*(^* - ^'*^^^^)■ (7) 

The mass eigenvalues are given by 



M 

*1,2 



1 



2 
with the mixing angle 



iML+MU T J{Ml^ +Ml) + A\Ml\ (8) 



tLL tRR' ~ y ^ tLL taR 



-M? 
cosOi = *^^ (9) 

'\M? \^ + (M? -M? '" ^ ' 

' tLR I ^ tl tj 



^2 
^tLL' 



3 



Mf -M? 
sinOi = *" '' (10) 

In the following we will consider a particular choice of the MSSM parameters that are 
defined, in the notations of PYTHIA6.4, in the following way: 

Mq = 270 GeV; M^ = 270 GeV; At = -500 GeV (stop trilinear coupling); 
tan(3 = 5; fi = -370 GeV; Mi = 80 GeV; M2 = 160 GeV. 

Note that in PYTHIA6.4 M^ corresponds to Mj (left squark mass for the third 
generation) and M^ corresponds to Mj . These parameters give Mj^ = 167.9 GeV, 
M + = 159.2 GeV and M^o = 80.9 GeV. This parameter point is compatible with all 
experimental data. We have chosen this value of Mj rather close to the mass of the 
top quark Mtop = 170.9 ± 1.8 GeV [15]. This means that one expects a rather large 
contribution from the top background, therefore, this choice makes the analysis most 
difficult. Finding a suitable set of cuts separating stop and top events is crucial. 

In general, the cross section for stop pair production at a fixed energy depends on the 
mass of the stop quark and the mixing angle 6^. Since the couplings of the Z^ to the 
left and right components of the stop are different, the cross sections depend significantly 
on the beam polarizations (see |7j, [8], [22]). By choosing appropriate longitudinal beam 
polarizations it is possible to enhance the cross sections. For example, for an electron beam 
with 90% left polarization the cross section would be larger than the unpolarized cross 
section by approximately 40%, for cos Oi = - 0.81 corresponding to the parameters given 
above. If in addition the positron beam has 60% right polarization, then the cross section 
is enhanced by approximately a factor of 2 compared to the unpolarized cross section. 
We note that a rather precise determination of the stop mixing angle 9i is possible by 
measuring the left-right asymmetry. The cross section for top pair production has also a 
characteristic dependence on the beam polarizations [22]. For example, the polarization 
of both beams leads to an increase of the cross section by about a factor of 1.5. 

3 Distributions of kinematical variables in stop and 
top production. 

In the following we present some distributions of different physical variables based on 
the sample of 5 ■ 10^ stop pair production events generated by PYTHIA6.4 and CIRCEl 
weighted with the electron-positron luminosity. Analogous plots are also given for 10^ 
generated background top events. In this Section all plots are obtained without any cuts. 

The ILC is a linear electron-positron collider with a center-of-mass energy of 200-500 
GeV and a high luminosity (peak luminosity of ~ 2 • 10^^cm~^s~^), upgraded to 1 TeV 
in the second phase. According to P] a total luminosity of 500 fb~^ is foreseen within 
the first four years of operation and 1000 fb~^ during the first phase of operation at 500 
GeV. A first run at i/i = 500 GeV will get a first measurement of the particle masses to 
optimize the threshold scan [5]. 

Fig. 2 a) demonstrates the total energy spectrum of the electron and positron beams 
Ele+ + E2e-, which is expected at 2Eheam = \/s = 500 GeV after taking into account 



beamstrahlung and other beam interaction effects (see, for instance [23]). Fig. 2 b) sliows 
tfie correlations of tlie beam energy fractions yi = Ei/Ebeam (i=e''', e~) of tlie colliding 
electron and positron beams. 



xlQ- 



-500- 



400 



300 



200 



100 



Entries 

Mean 

RMS 

Underflow 

Overflow 

Integral 



1000000 

494.3 

10.2 





1e+06 



°0 



a) 



100 200 300 400 500 600 700 800 9001000 

E1..+ E2„., GeV 



CM 1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0, 






Entries 
Mean x 
Mean y 
RMSx 
RMSy 
Integral 



1000000 

0.9758 

0.9758 

0.03829 

0.03807 

2.212e+05 



2498181280342 



221183248657 




j^ 



j^ 



j^ 



l1^ 



L±l. 



J^ 



J^ 



J^ 



J^ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



b) 



Vi 



Figure 2: a) total beam energy spectrum, b) beam, fractions correlations. 



3.1 Distributions in stop events. 

In Sections 3, 4 and 5 we will present plots for the kinematical distributions only for 
the energy -y/i = 500 GeV. Fig. 3 shows two kinematical distributions characteristic of the 
produced stop pair system, i.e., the number of expected events versus the stop transverse 
momentum PTj (plot a)) and its polar angle 6j (see plot b) (all in the e^e~ c.m.s.). As 
can be seen in Fig. 3 a) the stop transverse momentum PTj spectrum begins at PTj ^ 50 
GeV and has a peak near the kinematical limit. 




Mean 

RMS 

Underflow 

Overflow 

Integral 



158.2 

28.71 





2373 



50 100 150 200 250 300 350 400 450 500 

PT t„ GeV 



b) 











Mean 
RMS 


1.569 


£bU 






kPi 


0.5006 


50 
40 






Jl/^ 


Underflow 









/ ^ 


Overflow 
Integral 





-■- 


J 


J 


2373 


30 








/ 




20 


-■■ 


/ 










10 


r-- 


J 












0.5 


1 1.5 2 


2.5 . 


3 





Figure 3: Distributions of the number of events Nevent 
transverse momentum PTr ; b) stop polar angle Or ■ 
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a) stop 



In Fig. 3 and the following figures the vertical axis shows the number of stop and top 
events that may be expected for the integrated luminosity hint = 1000 fh^^. Taking 
the integral of the distributions one can get the total number of events expected for this 
luminosity. These numbers are shown as "Integral" values in the Figures. 



To find the jets we use the subroutine PYCLUS of PYTHIA with the distance measure 
utilized in the "Durham algorithm". The parameters of this jet finder, which is widely 
used in e+e" physics, are chosen such that the number of jets is exactly four (see also 



3.2 Distributions of the jets from W decay. 

According to the decay chain (2), the final state has to contain two jets due to the 
decay of one W boson into two quarks W ^ qi + qj (see Fig.l). 
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Figure 4: Energy spectra of the quarks from W boson decay, a) stop pair production; b) top 
pair production. 

Plot a) (for stop pair production) and plot b) (for top pair production) of Fig. 4 show 
the distributions of the energy Ew-quark of the quarks produced in the W boson decay 
(which we call " W^-quarks" ) . The first spectrum begins at zero and goes up to 140 GeV, 
with a mean value of 52 GeV, while the second spectrum begins around 8 GeV and goes 
up to approximately 200 GeV, with a mean value of 72 GeV. 
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Figure 5: PT spectra of the quarks produced in the W boson decay, a) stop pair production; b) 
top pair production. 



Figure 5 shows the transverse momentum PTw-quark spectra of the quarks produced 
in the W boson decay for stop (plot a)) and top (plot b)) production. The shapes of the 



PTw-quark spcctia of these "VT-quarks" are rather similar to the Ew-quark spectra. In 
the case of top production the " ly-quarks" are shghtly more energetic and have a larger 
transverse momentum than those from stop pair production. 

As the next step we take into account the hadronization of the "PF-quark" into a jet 
which we call "jetw"- Figs. 6 and 7 show the energy Ej^ty^ and transverse momentum 
PTjetw distributions of the corresponding "ly-jets". Plots a) and b) are for stop and top 
production, respectively. According to our choice of PYCLUS jet finder parameters there 
are two "jet^y" in the event. 
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Figure 6: Ej^ty^ energy spectra, a) stop pair production; b) top pair production. 



Comparing plot a) of Fig. 6 for the energy distribution of " W-jets" in stop production 
with the plot a) of Fig. 4 one observes that the corresponding mean value of the "W-jets" 
energy Ej^tw* ("^^ ^^e the notation W* in the stop case, for some details see below) in 
Fig.6 is about 16 GeV lower than the mean energy Ey/_quark of " VT-quarks" . It is also seen 
(plot a) of Fig. 4) that the peak position of " H^-quark" energy distribution {E^°_ 



quark 



25 



GeV) is shifted to the left by about 17 GeV {E^^^^^ ^ 8 GeV) when passing to the jet level 
(see plot a) of Fig.6). The end point of the Ejetw* distribution in stop case is somewhat 
lower than that one for the corresponding quarks. 
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Figure 7: "jetw " PT- spectra, a) stop pair production; b) top pair production. 

Analogously, the mean value and the peak position of the distribution of the transverse 
momentum of the " Vr*-quarks" , PTw-quarki shown in Fig.5 a), decrease by about 13-15 



GeV when passing to the jet level (see Fig. 7 a)), while the end point of the PTjetyy* 
distribution is a bit lower than the end point of PTw -quark distribution. 

Due to the different kinematics in top production mentioned above, the energy Ej^tw ' 
its peak position and the mean value of the " jetvK" energy distribution in the top case are 
practically equivalent to the Ey/_quark spectrum, peak position and the mean value of the 
corresponding "Pi^-quark" energy distribution (see Fig. 4 b) and Fig. 6 b)). Analogously, 
by comparing plots b) of Figs. 5 and 7 for PTw -quark and PTjetyy, one can see that the 
transverse momentum distribution in top production is stable under hadronization. 
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Figure 8: The invariant mass of two quarks M\y = Minv {guar kl + guar k2) , reconstructed from 
the vectorial sum of 4-fnomenta of two quarks that are produced in W ^ gi + qj decay, a) stop 
pair production; b) top pair production. 



Figure 8 shows the spectrum of the invariant mass Minv{quarki + quark j 



Ml 



w 



reconstructed from the vectorial sum of 4-momenta of the two quarks produced in W 
decay W — )■ g^gj. Plot a) is for stop pair production, plot b) is for top pair production. 
In plot a) of Fig. 8 one clearly sees the virtual nature of the W boson in the stop pair 
production case which we denote by a (*) in W* . Hence, in the stop case the invariant 
mass of two quarks produced in the decay of the virtual W* is smaller than the mass of 
a real W boson. In top production (see plot b) of Fig.8) there is a peak in the invariant 
mass distribution at the mass value of the real W boson. 
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Figure 9: Number of generated events versus the reconstructed invariant mass of "all-non-b- 
jets". a) stop pair production; h) top pair production. 



Figure 9 shows the corresponding plots at the jet leveL The invariant mass in Fig. 9 
is built of two "jetw" - jets produced in the W boson decay (or, shortly, '"JETSw")- 
One can see from plot a) that in the stop case the spectrum of the invariant masses 
M]y* = Minv{JETSw) is strongly shifted to the left. It has the main peak in the 
interval from 20 to 30 GeV and two tails. The left tail is rather short while the right one 
is very long and it spreads up to Mw* ~ 180 GeV. There is another not so high peak of, 
approximately, 65% of the height of the main peak, which is seen near the position of the 
peak shown in the quark level plot a) of Fig. 8 near the point Mw* = 80 GeV. As seen 
from plot b), in the top case the position of the W-peak at the jet level remains (with a 
high precision) at the same value of Myy as at the quark level in plot b) of Fig.8. There 
appear a shift in the mean value (about 6 GeV) and some tails on both sides. 

3.3 6-quark and 6-jet distributions in stop and top production. 

In the case of stop decay into a 6-quark and a chargino, ti — )■ bxf, the jets produced 
in 6-quark hadronization are observable objects. Their features are interesting from the 
viewpoint of experimentally distinguishing the stop signal events from the top background. 
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Figure 10: b- and b-quark energy spectra, a) stop pair production; b) top pair production. 

In Fig. 10 we show in plot a) for stop pair production and plot b) for top pair produc- 
tion the distributions of the energies Eh of the b- and 6-quarks (which we do not distinguish 
in the following) produced in stop and top decay chains (2) and (3), respectively. Both 
spectra begin at £'b ~ 4 GeV, corresponding to the b-quark mass, but look very different. 
The b-quark energy spectrum in stop production strongly increases up at ii^h ~ 4 GeV and 
strongly decreases at Ef, ^ 17-19 GeV. The corresponding spectrum in top production is 
much harder and its main part is concentrated within the interval 25 < Ei, < 170 GeV. 
The mean values of the b-quark energies are about 11 GeV and 94 GeV in stop and top 
production, respectively. This means that in the stop case (see plot a)) the b-quark takes 
a smaller part of the stop energy Ej ^ 250 GeV than the b-quark gets in the background 
top case (see plot b)). 

Figure 11 shows the transverse momentum PTf, spectra of 6-quarks for stop (plot a)) 
and top (plot b)) production. Comparing plot a) of Fig.ll with plot a) in Fig.3, one 
can conclude that in stop pair production the 6-quarks have only a small fraction of the 
transverse momentum of the parent stops. The shape of the PTf, spectrum of 6-quarks in 



9 



the stop case (see plot a) in Fig. 11) is similar to the shape of the Eb spectrum (see plot 
a) in Fig. 10). This means that in the stop decay the transverse component of the 6-quark 
momentum is larger than the longitudinal component. 
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Figure 11: b- and b-quark PT spectra, a) stop pair production; b) top pair production. 

The kinematical distributions of the 6-quarks in top decay are quite different. As seen 
from plots b) of Fig. 10 and Fig. 11, the 6-quarks produced in the top decays are very 
energetic. Most of the top events have E^ > 25 GeV and PTf, > 20 GeV. The difference 
to stop decay is easily understandable. The stop decays into a heavy chargino, whereas 
the top decays into a real W boson whose mass is only half of the mass of the chargino 
M ±. Therefore, the 6-quarks in top decays have a larger phase space than the 6-quarks 
m stop decays. 
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Figure 12: a) b- and b-quark polar angle @b distribution in stop production; h)cos{b,b) distri- 
bution in stop production. 

The distribution of the polar angle B^ of the 6-quarks in stop production is presented in 
plot a) of Fig. 12. Plot b) of Fig. 12 contains the cos{b, h) distribution, where cos{b, b) is the 
cosine of the opening angle between the 3-momenta of the b- and 6-quarks produced in the 
same stop event. It demonstrates that most of the b- and 6-quarks move in approximately 
opposite directions, but some are in the same hemisphere. Thus in the experiment we 
may expect a similar angular distributions of the corresponding b- and b- jets. 
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As the next step, we take into account 6-quark hadronization into a 6-jet. Technically, 
6-jets are defined as jets that contain at least one B-hadron. Their decay may be identified 
by the presence of a secondary vertex 



Figs. 13 and 14 show the energy -Eb_jet and transverse momentum PTh-jet distributions 
of the corresponding 6-jets. Plots a) and b) are for stop and top production, respectively. 
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Figure 13: b-jet energy spectra, a) stop pair production; b) top pair production. 
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Figure 14: b-jet PT- spectra, a) stop pair production; b) top pair production. 

Comparing plots a) of Figs. 13 and 14 for the b- and 6-jet energy Eb-jet and jet trans- 
verse momentum PT^-jet distributions in stop production with the Figs. 10 a) and 11 a), 
one observes that the corresponding mean values of jet energy and transverse momentum 
are about 15-17 GeV higher than those for the quarks. The reason for this is that in the 
stop case the end points of the energy distributions for the 6-jets and 6-jets are higher 
than those for the corresponding quarks due to appearance of long tails at higher E^- jet- 
It is worth emphasizing that at the same time the position of the peak in 6-jet energy 
distribution Eb-jet (see plot a) of Fig. 13) in the stop case is only slightly shifted to higher 
values (by about 3-4 GeV) compared to the Ei, peak in plot a) of Fig 10. Also note that 
there are practically no changes of the peak positions of the quark level PT^ and the jet 
level PTb-jet distributions shown in plots a) of Figs. 11 and 14, respectively. Thus we can 
say that the peak positions of energy and transverse momentum distributions turn out to 
be rather stable when passing from 6-quark to 6-jet level. 



11 



Now let us turn to the analogous energy and transverse momentum distributions 
obtained for the case of top background. As seen from comparison of the energy plots 
b) of Figs. 10 and 13, as well as the transverse momentum plots b) of Figs. 11 and 14, 
the mean values of the 6-jet and 6-jet energy and transverse momentum distributions are 
about 1-3 GeV smaller than the mean values of the corresponding 6-quark and 6-quark 
distributions. 

Let us summarize the results which were obtained in subsections 3.2 and 3.3 by the use 
of PYCLUS jet finder. First, it was found that in the case of top background production 
the characteristic parameters of energy and transverse momentum distributions of jets 
stemming from W decay and of 6-jets, produced in 6-quark hadronization, practically do 
not differ from the parameters of their parent quarks distributions. 

This picture changes quite noticeably when we consider the case of stop production 
with its further decay through the channel ti — )■ 6xi • I^ this case the 6-quarks are much 
less energetic than the b-quarks produced in top decay t — t- hW^. It was observed that 
the use of the same PYCLUS jet finder in the stop case leads to a noticeable redistribution 
of jet energies and, correspondingly, of their transverse momentum. Namely, the mean 
values of jet energy Ej^t^^ and jet transverse momentum PTjet^, are about 15-17 GeV 
smaller than the energy -EvK-guarfc and transverse momentum PTy/_quark of parent " W- 
quarks" (stemming from W boson decay), while the mean values of 6-jet energy Eh_jet 
and jet transverse momentum PTf,_jet are by about 15-17 GeV higher than the energy Eh 
and PTh of parent 6-quarks. 

In the following we shall return to this subject and consider the set of physical vari- 
ables which shall take into account this effect of energy redistribution in a case of stop 
production. 

4 Distributions of the signal muons. 



To select the signal stop pair production events shown in the left diagram of Fig.l 



one has to identify the muon from the W decay. The corresponding energy E, 
transverse momentum PT, 



sig—mu 



and 



sig—mu 



distributions of the signal muons are shown in Fig. 15. 
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Figure 15: a) Energy distributions of signal muons. b) PT distributions of signal muons. 
There are, however, also muons in the event coming from leptonic and semileptonic 
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decays of hadrons. Fig. 16 a) and b) show, respectively, the energy Edec-mu and the trans- 
verse momentum PTdec-mu spectra of these background muons stemming from hadron 
decays within the detector volume (for which we took the size from |3], [5]). It can be 
seen that the decay muons have a rather small energy Edec-mu and transverse momentum 
PTdec-mu- Their mean values are about 0.79 and 0.63 GeV, respectively. The analogous 
spectra for the signal muons in Fig. 15 show that the signal muons have a much higher 
energy Egig-mu and transverse momentum PTsig-mu- The mean value of the signal muons 
energy -E"™^"^^ = 47.6 GeV is about 60 times higher than the mean value of the energy 
of the decay muons. An analogous difference can be seen between the mean values of 
transverse momenta PT of signal and decay muons. One can cut off most low-energy 
decay muons rejecting those with Emu < 4 GeV. Such a cut leads to a loss of about 15-20 
signal events as seen from the plot a) of Fig. 15 (the bin width in this plot is 2 GeV). 
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Figure 16: Distributions of background muons. a) Energy distribution; b) PT distribution. 

We have also studied another way to select the signal muon from W decay. If the axes 
of all four jets in the event are known, then in general the signal muon has the largest 
transverse momentum with respect to any of these jet axes. 

5 Global variables useful for background separation. 

In stop pair production the two neutralinos and the energetic neutrino from the 
W boson decay escape detection. The simulation with PYTIIIA6 allows us to estimate 
the missing energy and the missing transverse momenta that are carried away by these 
particles. We also take into account the non-instrumented region around the beam pipe 
given by the polar angle intervals < 7° and > 173°. 

The distributions of the total missing Emiss-tot energy for stop production and top 
production are presented in a) and b) plots of Fig. 17, respectively. In stop pair produc- 
tion, see plot a), the Emiss-tot spectrum starts at 200 GeV. In top pair production (plot 
b)), where two neutralinos are not present, the missing energy Emiss-tot is much smaller, 
going from ^ 10 GeV to ^ 260 GeV. 

Figure 18 shows the distributions of the total visible energy in the event, E^s-tot^ 
in stop production (plot a)) and in top production (plot b)). The large missing energy 
in stop production (Fig. 17) is related to the low visible energy (Fig. 18), while in top 
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Figure 17: Missing energy Ejyiiss-tot distribution, a) stop pair production; b) top pair produc- 
tion. 

production the low missing energy correlates with the large visible energy. A cut on 
the total visible energy of approximately Ems-tot < 220 GeV g would eliminate almost 
completely the top background, while leaving the most part of the signal events. 
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Figure 18: Total visible energy Ems-tot distribution, a) stop pair production; b) top pair 
production. 

Another useful observable is the scalar sum of the moduli of the transverse momenta in 
an event PTgcaisum = Y.i^=r 1-^^*1) where the sum goes over all (A^p"^*) detectable particles 
in the event. Figure 19 shows the distributions of the scalar sum of the transverse momenta 
for the stop production (plot a)) and for top production (plot b)). It is seen that the 
restriction PTscaisum < 150 GeV would lead to a good separation of the stop signal events 
from the top background. 

Let us consider also the invariant mass of the system that contains all observable 
objects in the final state. This invariant mass is the modulus of the vectorial sum of the 
4-momenta P"^ g of all (A^-'^* = 4, n = 1, 2, 3, 4) jets in an event plus the 4-momentum of 



the signal muon P. 



A" 



^That is equivalent to setting a lower limit for the missing energy. 
^ The 4-momentum P" ~ [Pq ^ P")i of the n-th jet includes its energy Pq 
of the 3-dimensional momentum P" = {P^,P^,P^). 
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Figure 19: PTg^alsum- distribution, a) stop pair production; b) top pair production. 



Minv {AUjets, mu) 






(11) 



The distribution of this invariant mass is shown in Fig. 20. Plot a) shows the results 
for stop pair production while plot b) is for top pair production. As seen from these plots, 
the cut Minv { All jets,fi) < 200 GeV will give a good separation of signal stop and top 
background events. 
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Figure 20: Distribution of number of events, Nevent {Lint = 1000/6~^), versus the reconstructed 
invariant mass of all jets and signal muon Minv{Alljets,mu). a) stop pair production; b) top 
pair production. 



6 Global jet variables and cuts. 

An even more efficient separation of the signal and the background can be obtained 
by using the invariant mass of all jets Minv{Alljets) which is the modulus of the vectorial 
sum of the 4- moment a of all N^"^^ jets in an event 



Miny {All jets) 



^^ pn \2 



(12) 
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The distribution of this invariant mass is shown in Fig.21. Plot a) presents the result 
for stop pair production while plot b) that for top pair production. It is seen that the 
application of the cut Min^^AUjets) < 160 GeV leads to practically a complete separation 
of signal stop and top background events. 
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Figure 21: Distribution of number of events, Ng^ {Lint = WOO fb ^), versus the reconstructed 
invariant mass of all jets Minv{Alljets). a) stop pair production; b) top pair production. 

Another variable that can also be used for the separation of the signal and the back- 
ground is the "missing" mass Mmissing (we use a/s = 500 GeV) 



M 



missing 



, AT jet , AT jet 

(V^ - (E.., E^t + E,)? - (E„.i Pie. + P 



(13) 



This variable takes into account the contribution of those particles that cannot be 
registered in the detector (neutrinos and neutralinos) . The distributions of this invariant 
"missing" mass are given in Fig. 22. Plot a) shows the results for stop pair production while 
plot b) is for top pair production. As can be seen from these plots, the cut Mmissing > 250 
GeV also allows us to get rid of most of the top background contribution. 
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Figure 22: Distribution of number of events, N^v {Lint = 1000/6 ^), versus the missing mass 
variable, a) stop pair production; b) top pair production. 



16 



6.1 Cuts and signal-to-background ratio. 

To diminish the influence of jet energy redistribution effect, discussed in subsections 

3.2 and 3.3, we shall use the cuts considered above for Mmissing and Minv{AUjets). These 
variables, by definition, include the total 4-niomentuni of all jets, defined as the vectorial 
sum of the 4-momenta of all jets. Therefore they do not suffer from energy redistribution 
between jets. Consequently, we use the following three cuts to separate the signal and 
background events: 

• there must be at least two 6-jets in an event: 

Nb-jets > 2; (14) 

• the invariant missing mass must be larger than 250 GeV: 

Mmissing > 250 GeV; (15) 

• the invariant mass of all jets must be smaller than 160 GeV: 

Minv{Alljets) < 160 GeV. (16) 

These three cuts reduce the number of top background events from ^ 3.5 * 10^ to 
Ri 12 and the number of stop signal events from 2373 to 1806. So, the cuts improve the 
signal-to-background ratio from S/B = 0.066 to S/B ^ 143 losing about 23.8% of the 
signal stop events. The efficiency values for the cuts (14)-(16) are calculated for y/s = 500 
GeV. We define them as the summary efficiencies. It means that if ei is the efficiency of 
the first cut (14), £12 is the efficiency of applying the first cut (14) and then second cut 
(15). Analogously, £123 is the efficiency of the successive application of the cuts (14), (15) 
and (16). The following results are obtained: 

For SIGNAL STOP events: €1 = 0.84; £12 = 0.78; £123 = 0.76; 

For BACKGROUND TOP events : ei = 0.94; £12 = 0.001; £123 = 3.3 ■ 10"^ 



6.2 Signal and background cross sections. 

We give the cross sections and the numbers of events for stop pair production and 
top pair production for five energies y/s in Table 1 without cuts and Table 2 with cuts. 
To calculate the rates of stop production at different energies we take the value of the 
integrated luminosity Lint = 1000 fb~^ for all energies. 

It is seen that the cuts (14)-(16) lead to a strong suppression (more than 3 orders of 
magnitude) of the background top contribution and a moderate loss of signal stop events 
produced in the energy range 400 < y/s < 800 GeV. 

For the stop mass chosen, the largest number of signal events is expected at a/s = 500 
GeV. Let us note that, according to [21], a 50% efficiency of the separation of b -jets 
and 80% of the corresponding purity can be expected. It means that in reality to get 
1806 signal stop events, reconstructed with an account of b- tagging, we will need about 
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2.5 times higher statistics than that provided by the integrated luminosity 1000 fb ^ at 
^ = 500 GeV. 

It is worth noting that, as seen from Table 2, even with the use of the fixed parameters 
of the cuts (14)-(16) the number of signal stop events that can pass these cuts grows 
rapidly with the energy in the region 400 < ^/s < 500 GeV where mass measurements 
will be done in the first phase. These measurements may allow to enlarge the sample 
of collected signal events and to perform a precise measurement of the stop mass. The 
region 500 < ^/s < 800 GeV, which will be available in the second phase of ILC operation, 
allows to gain a much higher stop statistics, as seen in Table 2. A complete analysis based 
on adjusting the parameters of the selection cuts (15)-(16) for each of the above energies 
intervals will be presented in forthcoming papers. 

Table 1: The cross sections and the number of events for stop and top pair production 
before cuts. Lint = 1000/6"^ is assumed for each energy 
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Table 2: The same as in Table 1 but after applying cuts (14)-(16). 
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7 Determination of scalar top quark mass. 

Another variable of interest is the invariant mass Minv{b — jet, JETSw)'- u 



M,n.ib-jet,JETSw) = JiPt 



b-jet 



PjETSwY-: 



(17) 



which is constructed as the modulus of the vectorial sum of the 4-momentum Ph-jet of the 
6-jet, plus the total 4-momentum of JETSw system, i.e., non-6-jets stemming from the 
W decay {Pjetsw ~ Pjetiw + Pjet2w^ ^^ there are only two jets allowed to be produced 
in W decay). More precisely, if the signal event contains a ix~ as the signal muon (see 
Fig.l), we have to take the 6-jet (6-jet in the case of /z^ as the signal muon). This is 



We follow here the notations of subsections 3.2 and 3.3 



only possible if one can discriminate between the b- and 6-jets experimentally. Methods 
of experimental determination of the charge of the 6-jet (&-jet) were developed in |21] . 
In the present paper we do not use any b-tagging procedure. The PYTHIA information 
about quark flavor is taken for choosing the b- and 6-jets. 

The distributions of the invariant mass of the "b-iet+JETSw" system is shown in 
Fig.23. Their analogs Minv{b, 2 quarksw), obtained at quark level, are presented in Fig.24. 
Plots a) of these two Figures show the results for stop pair production while plots b) are 
for top pair production. These distributions are obtained without any cuts. 
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Figure 23: The spectra of the invariant masses Minv(h 
(16). a) stop pair production; b) top pair production. 
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Figure 24: The spectra of the invariant masses M^, 
a) stop pair production; b) top pair production. 
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In the top case the invariant mass Miny{b,2 quarksw) of the system composed of a 
6-quark and two quarks from W decay should reproduce the mass of their parent top 
quark (see Fig.l). The distributions of events dN^'"'^^^ jdMin^l^ GeV expected in each bin 
of 5 GeV versus the invariant mass Minv{b, 2 quarksw) of the parent three quarks as well 
as the invariant mass of jets produced by these quarks, i.e. Minvib — jet, JETSw), are 
shown for jet and quark levels in plots b) of Fig.23 and 24, respectively, for an integrated 
luminosity of 1000 fb~^. These distributions have an important common feature. Namely, 
they show that the peak positions at quark level and at jet level, practically coincide 
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to a good accuracy with each other as well as with the input value of the top quark 
mass Mtop = 170.9(±1.8) GeV. It is also seen from plot b) of Fig.23 that the quark 
hadronization into jets leads to a broadening of very small tails which are seen in the 
invariant mass distribution at quark level (plot b) of Fig.24). The right tail, which 
appeared at jet level (see plot b)) of Fig.23) is a bit lower than the left tail and is longer 
than the left one. One may say that the peak shape at jet level still looks more or 
less symmetric. The main message from these plots is that the appearance of tails due 
to quark fragmentation into jets does not change the position of the distribution peak, 
which allows us to reconstruct the input top mass both at quark and jet levels. 

An analogous stability of the peak position at the jet and quark levels for the stop case 
can be seen in the plots a) of Figs. 23 and 24. Note that, according to the stop decay chain 
(2), the right edge of the peak of the invariant mass distribution of the "6 + 2quarksw" 
system (see plot a) of Fig.24) corresponds to the mass difference Mj — M^o. 

Now let us recall that according to subsection 6.1 the application of the cuts (14)- 
(16) leaves only 12 background top events and saves about 76% of signal stop events. It 
means that the picture shown in plot b) of Fig.23 would change drastically and resemble 
a random distribution of the twelve top events in a rather wide interval U- In the case 
of stop pair production the distributions of the invariant mass of the "6 + 2quarksw" at 
quark level and of the system "b-]et+JETSw" at jet level system are shown in Fig. 25 
taking only those stop events that have passed the cuts (14)-(16). 
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Figure 25: The spectra of the stop signal events ({Lint = 1000/5 ^)) after the cuts (I4) - (16), 
versus the invariant mass Minv (b-jet, JETSw*)- a) at quark level, b) at jet level. 



It is seen that the peak position of the stop distribution at jet level Mj„^(6-jet, 
JETSw*), obtained after application of the cuts (14)-(16) (plot b) of Fig.25), coincides 
with the peak position at quark level (plot a) of Fig.25) as well as with the peak positions 
in plots a) of Figs. 23 and 24 obtained without any cuts. Let us note that the observed 
stability of the peak position in both plots of Fig.25 is due to the rather moderate loss 
of the number of events in the peak region (they change from ?» 550 to ~ 290) while 
passing from quark level to jet level. The interval 150-350 GeV in the plot b) of Fig.25 
can be used to calculate the width between the grid dots in this plot. It is found to be 
about 7.4 GeV. This number allows to estimate the position of the right edge of the peak 



^ for more details see Section 8 
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of Mi„^,(6-jet, JETSw) distribution, which seems to be shifted to the left side from 100 
GeV point by the distance which is a bit less than two dot intervals, i.e., by a bit less 
than 14.8 GeV. Thus, we can estimate that the right edge of the Mj„^,(6-jet, JETSw) 
distribution peak lies a bit higher than 85.2 GeV. 

Some remarks about the tails in the stop distributions are in order now. The origin of 
the right and left tails of the distribution shown in the plot a) of Fig. 25 can be clarified 
by the results of the stop mass reconstruction by calculating its invariant mass at quark 
level Minvip, 2 g-uar/csiy*, Xi) as the modulus of the sum of the 4- momenta of all three 
quarks and the neutralino (see Fig.l) in stop decay. These results are given in plot a) of 
Fig. 26 which shows a very precise reconstruction of the input stop mass at quark level 
withing the 5 GeV width of the bin containing the peak. Comparing plot a) of Fig. 25 
with plot a) of Fig. 26 one can conclude that at quark level the long left tail as well as the 
very small right tail in the distribution of M(„^(6, 2 quarksw) disappear when neutralino 
4-momentum is added to the 4-momentum of the " 6 + 2quarksw" system. 
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Figure 26: The spectra of the stop signal events after cuts ({Lint = 1000/6 ^)) versus the 
invariant mass Minv (b-jet, JETSw, Xi)- ^) ^^ quark level b) at jet level. 

The influence of the effect of the hadronization of the 6-quarks and of the quarks 
from W decay into jets is shown in plot b) of Fig. 26. This plot demonstrates that the 
hadronization of quarks into jets does not change the position of the stop mass peak, 
which practically coincides with the input value Mj = 167.9 GeV. It is also seen that the 
hadronization results in the appearance of more or less symmetrical and rather suppressed 
short tails around the peak position. The mean value of the distribution in plot b) of 
Fig. 26 is also very slightly different from the mean value of the quark level distribution 
shown in plot a) of Fig. 26. Both of these mean values are in a good agreement with the 
input value of the stop mass. It is worth mentioning that the shape of the peak in the 
stop plot b) of Fig. 26 looks similar to the shape of the peak in the top plot b) of Fig. 23 
which also demonstrates the stability (compaired to the quark level top plot b) of Fig. 24) 
of the reconstructed top mass peak position after taking into account the effect of quark 
fragmentation to hadrons. 

It is seen from plot b) of Fig.25 that by adding the mass of the neutralino M^o = 80.9 
GeV to the value of the right edge point of the peak Mj„i;(6-jet, JETSw) ~ 85.2 GeV 
one can get the lower limit for the reconstructed stop mass M- 
reproduces well the input value Mr = 167.9 GeV. 
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The simulation has shown (see plot b) of Fig.27) that the 12 background top events, 
which have passed our cuts (14)-(16), as discussed above, are mostly distributed in the 
region 30 < Mj„„(6-jet, JETSw)^ 140 GeV. This region is by more than twenty times 
wider than the 5 GeV width of the peak interval in the Minv{h-]et, JETSw) distribution 
which is shown in plot b) of Fig. 25 for the stop and which contains about 290 signal stop 
events left after the cuts. Therefore, we expect that in future measurements the contribu- 
tion of these twelve remaining top background events will not influence the position of the 
peak of the Mmt,(6-jet, JETSw) distribution (shown in plot b) of Fig.25) which allows 
one to reconstruct the input value of the stop mass by adding the mass of the neutralino. 

8 Results for top squark mass Mr = 200 GeV. 

In this section we want to discuss what will change if the mass of the top squark 
is different from the one we have chosen. In the present paper we have chosen a rather 
low scalar top quark mass (one of the lowest stop quark's masses that is allowed for the 
case of ti — > hxi decay channel). With increase of the stop mass the cross section for 
its production is decreasing. So, for example, for the case of Mj = 200 GeV and the 
integrated luminosity Li^t = 1000/6"^ the number of events per year at y/s =500 GeV 
is decreasing to 509 (after the cuts (14)-(16)). The mass Mj^ = 200 GeV is still below 
the highest allowed stop mass for the ti — )■ bxf decay channel (which is about 255 GeV) 
corresponding to M^+ = 159.2 GeV and M^o = 80.9 GeV. For stop masses below and 
above the described region, the stop will decay to other channels which we do not consider 
in this paper. The distribution of the invariant mass M(„^(6 — jet, JETSw) of the "6- 
iet+JETSw*" system for events which have passed the same cuts (14)-(16) is shown in 
Fig.27. Plot a) is for stop production with Mj = 200 GeV, plot b) is for top production. 
The top background also remains the same as it was given in Table 2 for y/s =500 GeV, 
i.e., about 12 events. 
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Figure 27: The spectra of the invariant masses Minv(b jet, JETSw) of the "b-jet+(aU- non-b- 
jets)" system after cuts (Lint = WOO fb^^). a) stop pair production; b) top pair production. 

The distribution in plot a) has a peak at Minv{b — jet, JETSw*) ~ HO GeV. One can 
also determine the mass of the stop quark following the procedure described in Section 7, 
but with less accuracy than in the case of the lower stop mass used in previous Sections. 
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9 Conclusion. 

We have studied stop pair production in electron-positron collisions within the frame- 
work of the MSSM for the total energies y/s = 350, 400, 500, 800, 1000 GeV. We assume 
that the stop quark decays dominantly into a chargino and a fe-quark, ii — )• bxf, and 
the chargino decays into a neutralino and a W boson, xf — )■ Xi^^^ where the W boson 
is virtual. One of the two Ws decays hadronically, W~^ — )■ qq, the other one decays 
leptonically, W~ — )■ /i~z/. 

We have performed a detailed study based on a Monte Carlo simulation with the 
program PYTHIA6.4 for -Eg+e- ~ V^ — 500 GeV and an integrated luminosity 1000 
fb~^. At this energy we expect the highest number of the signal events for the chosen 
stop mass M^^ = 167.9 GeV. The program CIRCEl is used to get the spectra of electron 
(positron) beams taking into account the effects of beamstrahlung. PYTHIA6.4 is used 
to simulate stop pair production and decay as well as top pair production being the main 
background. 

Three cuts (14)-(16) have been proposed to separate signal stop events and top back- 
ground events. For ^/s = 500 GeV and the luminosity 1000 fb~^ they give 1806 signal 
stop events with 12 background top events. This is different from the more complicated 
situation in stop pair production at LHC (see, for instance, [20]). 

We have shown that the determination of the right edge of the peak position of the 
invariant mass Minv{b-iet,JETSw*) distribution of the "6-jet+(all-non-6-jets)" system 
allows us to measure the mass of the stop quark with a good accuracy based on the 
statistics corresponding to an integrated luminosity of 1000 fb~^. For this the mass of Xi 
has to be known. 

As seen from the Table 2 the measurements at other energies in the regions 400 < 
i/i < 500 GeV and 500 < a/s < 800 GeV may allow to enlarge substantially the number 
of selected signal stop events and to perform a precise measurement of the mass of the 
scalar top quark. 

In the last Section we discussed the difference in the main invariant distribution for a 
Mj^ = 200 GeV. 

In conclusion we can say that the e^e~ channel together with the 77 channel considered 
in our previous Note [5] is well suited for the study of stop pair production at ILC. 
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